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ABSTRACT 

Recent studies found a correlation with ~3 (t significance between the local star formation measured 
by GALEX in Type la supernova (SN la) host galaxies and the distances or dispersions derived from 
these SNe. We search for these effects by using data from recent cosmological analyses to greatly 
increase the SNIa sample; we include 179 GALEX-imaged SNIa hosts with distances from the JLA 
and Pan-STARRS SN Ia cosmology samples and 157 GALEX-imaged SNIa hosts with distances from 
the Riess et al. (2011) Hq measurement. We find little evidence that SNe la in locally star-forming 
environments are fainter after light curve correction than SNe la in locally passive environments. We 
find a difference of 0.000±0.018 (stat-l-sys) mag for SNe fit with SALT2 and 0.029±0.027 (stat-|-sys) 
mag for SNe fit with MLGS2k2 (Ry = 2.5), which suggests that proposed changes to recent mea¬ 
surements of Hq and w are not significant and numerically smaller than the parameter measurement 
uncertainties. We measure systematic uncertainties of ~0.01-0.02 mag by performing several plausible 
variants of our an alysis. We find the g reat ly reduced significan ce of these distance modulus differ¬ 
ences compared to Rigault et al. (2013) and Rigault et al.|(2015l result from two im provements with 
fairly equal effects, our larger sample size and the use of jLA and 

5|). We find that 


criteria. Without these improvements, we recover the results of Rigault et al. 


both popula tions have more similar d ispersion in distance than found by Rigault et al. (2013), Rigault 
et al. (|2015|), and|Kelly et al.| (|2015|), with slightly smaller dispersio n tor locally passive (log[ lJsFR.) 
< —2 . 9 dex) SNe la ht with iVlLCS, the opposite of the effect seen by Rigault et al. (2015) and Kelly 


et al. (2015). We caution that measuring the local environments of SNe la in the future may require 
a higner-resolution instrument than GALEX and that SNIa sample selection has a significant effect 
on local star formation biases. 


1. INTRODUCTION 


Type la supernovae (SNe la) have been a key compo¬ 
nent in measuring the dark energy equatio n of state, w, 
with <6% uncertainty (Betoule et al.|2014) and the Hub- 
ble Gonstant, Hq, witn 3.37o uncertainty (Riess et al. 


2011 hereafter Rll). With such small error budgets, un¬ 
known systematic uncertainties affecting SNe la shape- 
and color-corrected absolute magnitudes could have seri¬ 
ous consequences for our understanding of dark energy, 
neutrino properties, and the global geometry of space. 

Although SNe la remain accurate distance indicators 
with ^10% uncertainty per SN, there are concerns about 
their ability to remain standardizable in galaxies that 
vary in mass, metallicity, star formation, age, and dust 
properties (e.g. S ullivan et al.||2010 Rigaul t et alJ[2MF 


Johansson et aE' 2013 " Childress et al. 2013r Even a 
small dependence of SN la luminosities on host galaxy 
properties may have a non-negligible effect on w due to 
the redshift evolution of galaxies or differences in sam¬ 
ple selection. Such an effect could also bias Hq due to 
the different demographics of Cepheid host galaxies com¬ 
pared to SNIa hosts. The lack of detection of such an 
effect at >3cr with samples of ^10^ SNe suggests that 


such effects are 5, 


10 % 

%/Ioo 


X 3 < 0.06 mag, or that they 


result from galaxy properties that are difficult to mea¬ 
sure robustly. These investigations are hampered by an 
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inability to define the nature of the SNIa correction a 
priori, complicating the interpretation of the significance 
of the correlations found a posteriori. If enough sources 
for a possible correlation are examined, a 3 ct result will 
always be found. 

The first widely accepted effect of host galaxy proper¬ 
ties on SNe la was confirmed by the detection of a ^0.07 
mag difference in mean corrected magnitude of SNe la 
with host masses >10^° Mf7). Identified by sev e ral inde¬ 
pendent stu dies includinglLampeitl et al. (2010), Sullivan 


et al. (12010 ), and Kelly et al. l 20lO), this effect has now 
been detected at >ba by |Hetoufe et af.| (|2014|) with a 
sample of 740 SNe la. 

Because it is unclear how the physics of a SNIa dis¬ 
tances could depend on its host galaxy mass, the most 
likely explanation is that host galaxy mass is merely 
tracing another physical property that could affect SN 
luminosity, such as m etallicity, stellar age, or dust. 


Dominguez et al. (20011 suggested that progenitor metal- 
licity could affect the SN luminosity by changing the 
Carbon-Oxygen ratio in the progenitor white dwarf, thus 
result ing in a lower Nickel m ass synthesized in the explo¬ 
sion. iHayden et al. (2013) found that a correction us¬ 


ing a star formation-based metallicity indicator reduced 
Hubble dia gram residuals more th an a simple host mass 
correction. Childress et al. (2013) found that dust and 
stellar age are also plausible explanations because they 
evolve with host galaxy mass. 

Different SN la progenitor ages could also exhibit sys¬ 
tematic differences in corrected magnitude due to the 
effects of metallicity or exp lo sion mechanism on ^^N i pro¬ 
duction (Maoz et al. 2014). Childress et al. (2014) sug- 
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gested that progenitor age could be the source of the 
host mass step, as older progenitors preferentially occur 
in non star-forming host galaxies. Because prog enitor 
age evolves with redshift, Childress et al. (2014) mod¬ 
eled a potential redshift-dependent bias in cosmological 
analyses. 

SNIa light curve fitters may also create biases by as¬ 
suming a universal relationship between color and ab¬ 
solute magnitude, independent of the dust composition 
of different SN la host s. Some preliminary e vidence has 
supported these ideas; |Scolnic et al.| (|2014b|) found that 
the correlation between blN ia color and absolute magni¬ 
tude has two different slopes for bluer and redder SNe, 
which may in part be due to dust properties. 

If the host mass step is indicative of one or more of 
these biases, galaxy properties in the vicinity of SN ex¬ 
plosions could be more strongly correlated with SN cor¬ 
rected magnitude than properties of the galaxies as a 
whole. Three recent studies used ~60—85 nearby SNe la 
to look at such properties and fo und that they affect 
the distances derived fro m SNe la. Rigault et al. (2013) 
and Rigault et al. (2015) found a correlation between lo- 
cal star formation and SN la Hubble residuals from the 
Nearby Superno va factory (Aldering e t al.|2002) and the 
CfA SN survey (Hicken et ai. 2U(J9a hereafter H09) by 
using the local star tormation rate density {T,sfr) to 
separate SNe la into those with locally passive (SN lae) 
and locally s tar-forming (SNIaa) environments. Rigault 
et al. (2015) (hereafter R15) found a mean difference m 
Hubble residuals between SNe lae and laa (hereafter re¬ 
ferred to as the LSF step) of ^0.09—0.17 mag at 2-4 (t 
significance with different light curve fitters. 

The fraction of SNe lae is different in the nearby 
Cepheid-calibrated SN la sample compared to the 
Hubble-ffow SNIa sample, and R15 found that SNe lae 
have mean corrected magnitudes ~0.15 mag brighter 
than SNe laa when fit with the MLCS light curve fit¬ 
ter and assuming the same Ry as the Rll Hq baseline 
analysis. They derived a correction to Hq: 


1 


= log(i7o) - -r)y- 


LSF bias correction 

( 1 ) 

where is the fraction of SNe lae in the Hubble- 

ffow SN sample and ip^ is the fraction of SNe lae in the 
Cepheid-calibrated sample. is the LSF step 

of 0.155 mag. By estimating (52.1±2.3%) and ip^ 
(7.0%), R15 estimate that the true value of Hq is reduced 
by ~3%. 

R15 also fou nd that SNe in highly star-forming r egions 


fit by MLCS dJha et al.||2007| 


Riess et al. 
lower dispersion in their Hubble residuals 


1996) have 
ran SNe in 


locally passive environments. |Kelly et al.| (|2015l came to 
the same conclusion by examining SiNela with high local 
star formation (Their Ssfr boundary is ~0.7 dex higher 
than the R15 lae/Iaa cut-off). R13 first found this effect 
using the SALT2 light curve fitter (|Guy et al.|2007|), but 
they could not r eproduce this result with H09 data. 

Both R15 and |Kelly et al.| (pM^ used GALEX FUV 
data to measure the star torrnation rate within a few kpc 
of SNe la positions. In this work, we use a similar method 
to examine whether the significance of the LSF step and 


reduced dispersion from SNe in locally star-forming host 
galaxies is reduced when we use the most current vintage 
SNe la distance estimates, use a much larger sample size, 
and vary the priors and assumptions used in the original 
analyses. 

Table shows th e sizes of the SN s ample s used in 
Rigaultetahl (j2013|), R15, |Kelly et al.| (|20i5|), and this 
work, along with the light curve htters used, the SALT2 
color param eters, and the MLCS prior on Ay- |Rigault| 
et al. (2013) used 82 SNfactory SNe with star formation 
estimated using loca l Ha from integral field spectroscopy. 
Rigault et al.|(2015D used ~100 SNe from the CfA3 sam- 
ple of HU9, w ith ^80 passing GALEX sample cuts, 
et al. (2015) used several surveys but made strict sample 
cuts and only used SNe with Hubble residuals < 0.3 mag, 
which would amount to a ^1.3cr cut for Rll data. 

By using a sample size ~2-3 times as large as those in 
the analyses above, we hope to obtain a robust measure¬ 
ment of the magnitude and uncertainty of the effect of 
local star-format ion on SNIa corrected magnitudes. §2 
presents our sample selection, and §3 discusses our LSF 
step and dispersion analysis. In §4 and §5 we present our 
results and discuss their significance, and our conclusions 
are in §6. 

2. DATA 

We used two samples of SNe for this analysis, one from 
the Rll measurement of Hq and the other from the dark 
energy equation of state measurements ofIBetoule et al. 
(I 2 M 4 I and Pan-STARRS (PSl;|Rest et al.|l2014t|Scolnic 
et al.[|2l)14a Scolnic et al. 2015, in prep). These two 


samples rely on many of the same SNe, but Rll use 
the MLCS lig ht curve fitter to p erform their baseline 
analysis while iBetoule et"^ (|2014|) a nd PSl use SALT2 
(jCuy et al.||2Ulu| |Betoule et al.||2014[ version 2.4). Each 
sampl e is ^2-3 times as large as the R15 and [Kelly et al.| 
(2015) GALEX-imaged host samples and removes the 
possibility of biases between our sample and the samples 
used in the most recent measurements of cosmological 
parameters. 


2.1. Riess et al. (2011) SNe 


The Hq determination of Rll use the MLGS2k2 light 
curve fitter for their baseline analysis. We use only their 
MLCS2k2 distance moduli, as JLA-l-PSl consists of a 
larger SALT2-fit SNIa sample with more robust light 
curve cuts and an updated SALT2 model and color pa¬ 
rameter, /3. The Rll sampl e consists of 1 40 SN e be- 
tween 0.023 < z < 0.1 from Hicken et al. (2009b) and 


Ganeshalingam et al. (2010). As one of the variants in 
their systematics section, Rll extend the lower bound of 
the redshift range to 0.01 after making peculiar velocity 
corrections (using res ults from Neill et al. 2007 and the 
Pike fc Hudson|2(^ dipole), giving 240 SNe (with pecu- 
liar velocity uncertainties added in quadrature to the dis¬ 
tances). Adopting this redshift range raises Hq by 0.8 km 
s“^ Mpc“^, or 0.26(7. We adopt this lower redshift limit 
of 0.01 as it allows us to add more SNe la to our sam¬ 
ple, although these nearby SNe have less weight in the 
likelihood approach outlined in §3 due to their included 
peculiar velocity uncertainties. In §4.1[ we examine the 
effect of restricting the redshift range to z > 0.023. Rll 
remove 4 (t Hubble diagram outliers but make no sample 
cuts based on light curve shape. Ay, or MLCS x^- 

























































































































3 



Figure 1. Hubble diagrams and Hubble residual diagrams for the JLA+PSl sample (SALT2 light curve fitter; left) and the Rll sample 
(MLCS light curve fitter with Ry = 2.5; right), with GALEX FUV-imaged hosts in red and SNe without GALEX FUV host images in 
blue. Out of a total of 249 SNe in the JLA+PSl sample, 207 were imaged by GALEX within 0.55 degrees of field center. In the Rll 
sample, 177 out of 239 SNe fit with Ry = 2.5 had GALEX FUV images. The MLCS data have slightly higher scatter, but both samples 
have intrinsic dispersions $0.2. 


Table 1 

Studies using local SF data 



SN Surveys 


SALT2 



MLGS 


SNe 

/^version 

0 

SNe 

/^version 

P(Ay) 

Ry 

Rigault+13 

SNfactory*^ 

82 

G07'= 

c 





Rigault+15 

CfA3 

77 

007*= 

2.48l°;l° 

84 

v0.06 

g-Av-/0.457 

1.7,2.5,3.1 

Kelly+15 

LOSS‘i,CfA2-4,CSP 




61 

v0.07® e 

-Av/0.3 = 0.02)^ 

1.8,3.1 

This Work 

CfAl-4,CSP,CTg,SDSS,SNLS,PSl 

187 

GIO*’ 

3.097 ± 0.062 

154 

v0.06 

g-Av-/0.457 

2.0,2.5,3.1 


Aldering et al.| )| 2002 [l. 


The value of ^ was blinded in |Rigault et aTj (j2013 . 

The Lick Observatory Supernova Search |Li~e^^L 20^^_;__^ 

® MLCS v0.07 used new spectral templates from |Hsiao et al.| | |2007[ |. This version was implemented in the SuperNova ANAlysis software jKessler e t al.| 
|2009b| SNANA). 

An exponentia l convolved with a n ormal distribution having <7 = 0.02 mag. 

^ Calan/Tololo ijHamuy et al.|l996[ |. 

et al.||201o| had improved uncertainty propagation and handling of residual scatter, a new SN la spectral energy distribution regularization scheme, 
and used a larger training sample with higher -2 SNe (see their Appendix A for details). 
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MLCS2k2 determines the distance modulus for each 
SN la by fitting for the light curve shape and extinction 
assuming an extinction prior and a value for the total-to- 
selective extinction ratio, Ry- Common extinction pri¬ 
ors include exponential distributions see Table 

[^, exponential distributions convolved with gaussians, 
a flat prior (with or without negative Ay allowed), and 
priors based on host galaxy information. Rll consider 
the latter two priors in their systematic uncertainty anal¬ 
ysis, and use an exponential with scale length 0.457 mag 
for their baseline analysis. Rll consider dust reddening 
laws with Ry = 1.5, 2.0, 2.5, and 3.1, using Ry = 2.5 
for their baseline analysis. Ry = 3.1 correspon ds to the 
Milky Way reddening law ( Cardelli et al |[T9M ). We ex¬ 
clude Ry = 1.5 from our analysis as such a low value is 
not t ypically used in cosmological analyses (e.g. Kessler 


tyi: ^_ 

et al.|200^ adopt Ry = 2.18±0.5 for SDSS cosmology); 


althou gh highly reddened SNe la tend to favor low values 
of Ry dBurns et al.|2014D, these SNe are usually excluded 
from samples used to measure cosmological parameters. 
H09, for example, use only SNe with Ay < 0.5. 

We queried GALEXj^for FUV images at the locations 
of these SNe, keeping only those with a angular distance 
from the held of view center (FOV radius) <0.55 deg to 
ensure accurate photometry and avoid rehection artifacts 
and distortion of the PSF near the detector edge. Of 
the 240 SNe used in Rll, we found 187 SN host images 
meeting this criterion, 157 of which remained after the 
sample cuts described in §3. A Hubble diagram of the 
Rll SN la sample is shown in Figure[l] There is less than 
0.01 mag difference in mean Hubble residual between the 
full sample and the GALEX-detected sample. No bias is 
expected for SNe with GALEX host images. 


2 . 2 . 


Betoule et al. (2014.) and Pan-STARRS SNe 


The most recent m easurements of w (Betoule et al. 
2014 Rest et aL]|2014 1 use the SALT2 light curve htter, 


and c ompute distance moduli using the equation (|Tripp 
19 ^: 


fi = m*g + axXi — /3xC — M, 


( 2 ) 


where ^ is the SN distance modulus, m*^ is the peak SN 
B band magnitude, Xi is the light curve stretch param¬ 
eter, and C is the light curve color parameter. SALT2 
adopts a linear relation between SN la color and luminos¬ 
ity with no prior. For consistency with the JLA cosmo¬ 
logical analysis, we only use the SALT2 fitter with these 
data. 

The nuisance parameters a, /3, and M (in this analy¬ 
sis, a single value independent of host galaxy mass) are 
simultaneously fit to the full supernova sample. In re¬ 
cent work, the value of (3 has risen due to changes in 
the SAL T2 model and l a rger S NIa samples. The value 
found by Betoule et al. (2014) is /3 = 3.102 ± 0.075, a 


difference of ^0.6 relative to the H09 value of 2.481° 
(used by R15). This could have an important impact 
on measuring the LSF step, which we discuss further in 
§5.1. In this analysis, we simultaneously ht JLA and PSl 
data together, findi ng /3 = 3.097 ± 0.062. In contrast to 
Betoule et al.|20l4|and following the R15 claim that the 


LSF step replaces the host mass step, we did not apply 
the host mass ste p in deriving this va lue. 

We limited the Betoule et al. (2014) Joint Light-curve 
Analysis (JLA) to z < 0.1 because the large GALEX 
PSF makes the star formation measurement non-local 
with FWHM ~8 kpc. This low-z sample i ncludes data 
from low-redshift surveys such as CfAl- 3 (Riess et al. 
19991 |Jha et al.||200 6[ iHicken et al.||2009bj), the Garnegie 
Supe rnova Project dHa muy et al.||2U06| IStr itzinger et al. 
2011|) and Calan/Tololo (IHamuv et al.|l996p, and surveys 

K. 


extending to higher z such as SUSS ( Kessler et al.|2009a 
25 SNe after sample cuts) and SNLS (|Gonley et al.|2011 
no S Ne after sample cut s). We added low-z GtA4 SJNe 
(|2012l used in the PSl analysis), PSl 


from Hicken et al. 


SNe from Rest et al. (12014 1 and the upcoming 4-year PSl 


cosmological analysis (12 SNe after sample cuts; Scolnic 
et al. 2015, in prep). For both JLA and PSl, p eculiar 
velocit ies are corrected follo wing Neill et al. (2007) based 


on the Hudson et al. (2004) modeT 


The cuts applied to these data are listed in [Betoule 
et al. (2014| their Table 6 and Appendix A). They make 
light curve shape, color, and SALT2 fit probability cuts 
(requiring a fit probability >0.01). We applied these 
same cuts to PSl SNe, and removed 3.5cr outliers from 
the full sample, inclu ding the 4 >3cr outliers removed by 
Betoule et al.| (2014|). 

The JLA and PSl samples with 0.01 < z < 0.1 con¬ 
tain a total of 249 SNe. 207 were found in GALEX with 
FOV radius <0.55 deg and 179 remained after the sam¬ 
ple cuts described in §3. We found no significant differ¬ 
ence (<0.01 mag) between mean Hubble residual of the 
GALEX-detected sample and the full sample. 

Figure[^shows a Hubble diagram for SNe in both sam¬ 
ples with and without GALEX imaging. Our cosmolog¬ 
ical fits used Om = 0.3, Ha = 0.7, w = —1, Hq = 70 km 
s“^ Mpc“^ and determined the absolute SN magnitude 
M from a least squares ht to the Hubble residuals. 

3. MEASURING THE STAR FORMATION DENSITY 

R15 used the following procedure to measure the local 
star formation density, Esfr, and its relation to SN dis¬ 
tance estimates. We summarize the principal st eps b e low 
and describe the differences in our analysis in )]3.1[ §3.2| 
discusses our systematic error treatment. Table[2[gives a 
summary of the quality cuts applied to our SN la sample 
and the number of SNe remaining after each cut. 

1. R15 measured GALEX FUV aperture photometry 
at the location of the SN using a 4 kpc aperture di- 
amet CT. They applied Mil ky Way dust corrections 
from Schlegel et al. (1998), where th e FUV extinc¬ 
tion A priv is 7.9 x h(B — V) (R15; |Cardelh et al. 
T989|. 


2. The photometry was corrected for host galaxy 
extinction in the FUV based on the measured 
FUV—NUV colors, which were converted to extinc¬ 


http://galex.stsci.edu/GalexView/ 


tion using the relation from Salim et al. (2007). 
A Bayesian prior of Apjjy = 2.0 ± 0.6 for star¬ 
forming galaxies was also applied (the hnal Apjjy 
was a weighted mean of the prior and the measured 
A-puv)- R-15 made no dust correction for passive 
galaxies. 

To determine whether each galaxy was globally 
star-forming or passive, they used Esfr measure- 
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Figure 2. Four host galaxies from our sample in SDSS gri images, with smoothed GALEX FUV contours marking the star-forming 
regions (log(SsFR) > —2.9) and the SNIa positions marked in red. Two apertures are overlaid, the local aperture size from R15 (4 kpc 
diameter) and the local aperture size from |Kelly et al.| (|2015| 10 kpc diameter). We assumed Apijy = 2.0 for the three star-forming 
galaxies. For the passive host of SN 2003ic, none ot the galaxy would be considered locally star-forming for Apiiv = 0? but we show dotted 
contours to indicate the effect of assuming 2 mags extinction. The 4 kpc diameter aperture appears to be a good approximation for the 
local star-forming environment while the 10 kpc aperture extends well beyond the local star formation environment for SN 2002ha and 
encompasses most of the galaxy for SN 2006en. Both the size of our local apertures and our prior on Apiiv have an important effect on 
our results, so we vary both in our systematic error analysis. 


Table 2 

SN Selection Cuts 


No. SNela 

JLA+PSl 

o 

oi 

II 

Rll 

=2.5 

II 

cii 

Initial Sample 

249 

240 

239 

237 

GALEX FUV data exist 

212 

189 

188 

187 

FOV radius < 0.55 deg 

207 

181 

180 

179 

Global SFR known 

207 

178 

177 

176 

Inclined SNe Removed 

179 

157 

156 

155 
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merits from Neill et al. (2009 Esfr > —10.5 is 


star-forming), who tit synthetic templates to the 
SN host UV-l-opt ical spectra l energ y distributions 


(SEDs). Because 
unavailable for 


Neill et al. (2009) SED fits were 
40% ot their hosts, R15 used mor¬ 


phology for these, treating galaxy types Sa and 
later as star-forming (a less accurate method). 


3. To minimize the effects of locally passive regions 
projected on top of locally star-forming regions (see 
R15, Appendix B.2), R15 removed SNe with host 
inclination angles >80° from their sample. 


4. Based on their photometric and dust correction 
uncertainties, R15 calculated the probability of a 
SNIa being above (P(Iaa)) or below (P(Iae)) the 
log(EsFR) = —2.9. 

5. R15 used a maximum likelihood approach (outlined 
in ( |A| ) to determine the difference in corrected mag¬ 
nitude and dispersion between SNelaa and lae. 


3.1. Our Analysis 

We largely used the same methodology as R15, but 
improved the following aspects of the analysis: 


1. We used the Schlafly fc Finkbeiner (2011) dust cor¬ 
rections instead ot the|Schlegel et al.|p998D correc¬ 
tions used by R15, resulting in a ~l47o reduction 
in our extinction valnes. 


2. We used SDSS NUV—r color instead of morphology 
as a diagnostic of global SFR when UV-boptical 
SED fits were unavailable. 


3. For SNe outside the isophotal radii of their host, we 
did not make a dust correction as we expect these 
SNe to be minimally affected by extinction. 

4. We made a slightly more conservative inclination 
cut, removing galaxies with inclinations >70°. 

5. Using our maximum likelihood model, we fit for 
both SN laa and SN lae dispersion when determin¬ 
ing the LSF step to allow for the possibility that 
these two quantities are significantly different and 
affect the magnitude of the step. 


We discuss our changes and methodology in further de¬ 
tail below. However, these ch ang es have only minor sig¬ 
nificance on our results (see (4.6). Our method of max¬ 
imum likelihood estimation for calculating the LSF step 
is described in detail in the Appendix. 


3.1.1. FUV Aperture Photometry 
We used the same baseline 4 kpc aperture diameter 


as R15 for our photometry but corrected for Milky Way 

FUV extinction using the 

Schlafly & Finkbeiner 

( 

2011) 

dust correction^ 

instead of the Schlegel et al. 

( 

1998) 

corrections usecT 

by R15. 

Schlafly V h'inkbeiner 

( 

2011) 

derive a ^14% correction tor the Schlegel et al. 

( 

1998) 

dust maps based 
SDSS stars. 

on the expected vs. measured col 

ors ot 


Using GALEX to estim ate local sta r form ation, as in 
Rigault et al.| (2015) and Kelly et al. (2015) is compli- 
cated by the large GALEX PST', 5.4" full width at half 



Figure 3. SDSS gri images of two spiral galaxies from our sample 
with SN positions marked in red and SExtractor-based isophotal 
radius estimates {R = 3) shown in blue. We corrected SN 1996bl 
for dust but did not correct SN 1999cw, as it exploded just outside 
the isophotal radius of its host galaxy and thus is beyond nearly 
all of its host galaxy’s dust. 

maximum (FWHM) in the NUV and 4.5" in the FUV, 
which serves as a lower limit to the size of the local re¬ 


gion that we can measure. Kelly et al. (2015) used a 10 
kpc a perture diameter to measure local star formation, 
while 
adopt 


Rigault et al. (2015) used a 4 kpc diameter. We 
the R15 4 kpc diameter in this work. 

Figure shows representative hosts from our sample 
with FUWbased log(EsFR) > —2.9 contours to demon¬ 
strate the size of these apertures relative to their star¬ 
forming regions. A 4 kpc aperture appears to be a rea¬ 
sonable approximation to the local SN la environment in 
these cases, while a 10 kpc aperture radius encompasses 
the majority of the SN 2006en host. In the case of SN 
2002ha, it is unclear whether either aperture is small 
enough to capture the star formation environment at the 
SN location. 


3.1.2. Host Galaxy Extinction Correction 

There are three principal differences between our local 
dust correction and that of R15. First , for galaxie s with- 
out star formation rates (SFRs) from|Neill et al. (2009) 
(45% of our sample), R15 used morphological iiitorma- 
tion to determine whether or not a galaxy was globally 
star-forming. However, GALEX NUV - SDSS r mag¬ 
nitude is a more reliable discri minator between p assive 
and star-forming galaxies (e.g. Salim et al. 2007 their 
Fig. 1). Passive galaxies have iNUV-r ^ 5, while star¬ 
forming galaxies have NUV-r < 4. For the 45% of our 
sample with SDSS images, we corrected for dust in galax¬ 
ies th at had NUV-r < 4.5 bas ed on SExtractor photom¬ 
etry (Berlin fc Arnou^|199 6 ). Fo r the final 19% of our 
sample without liN eill et al.l (|2U09|) SFR or SDSS images. 


we used morphology as an estiniiate of global star for¬ 
mation and performed a local dust correction for Sa and 
later-type galaxies. We removed 3 morphologically am¬ 
biguous hosts from our sample (SN 2005eu, SN 2006ah, 
and SN 2006is). 

Second, SNe la near the edges of galaxies should have 
negligible local dust. We used SDSS and, when n ecessary. 
Digitized Sk y Survey image^ to estimate the Sullivan 
et al.| (|2006|) SExtractor-based R parameter, which gives 
the SN separation from the host normalized by the size 
of the host galaxy. For the 28% of SNe approximately 
outside the isophota l radius of their host galaxy (R > 3; 
Sullivan et al. 2006), we did not correct for local dust 


® http://irsa.ipac.caltech.edu/applications/DUST/ 


http://archive.eso.org/dss/dss 
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regardless of the Salim et al. (2007) extinction estimate, 
which does not apply for passive, low-dust regions. R15 
dust-corrected all SNe in globally star-forming hosts, re¬ 
gardless of the location of the SN. Figure shows two 
examples of spiral host galaxies and their approximate 
isophotal radii. 

In total, our decision to apply or not to apply a dust 
correction was different from that of R15 for 14% of H09 
SNe (13/92 SNe). For 7 of these 13 SNe, we did not apply 
a dust correction because the SN was outside the isopho¬ 
tal radius of its host. The other 6 SNe had morphology- 
based SF classifications that disagreed with our NUV-r 
data. 

Finally, we adopted a slightly more conservative in¬ 
clination cut, removing galaxies wi th inclinations >70° 


based on the Tully & Fisher 
This removes an additional 


(19771 axial ratio method, 
b SlNe from the JLA-l-PSl 


sample and 11 from the Rll sample. In total, the incli¬ 
nation cut removes ~13% of our sample. 

3.2. Varying the Baseline Analysis 

For a robust result, we performed several plausible 
variants of our baseline analysis (R15 used a similar 
method to evaluate the robustness of the LSF step). We 
used the standard deviation of the measured LSF step 
from all variations to estimate our systematic error. 

Our FUV—NUV color measurements have a median 
signal-to-noise ratio of 3.02. Due to such large photo¬ 
metric uncertainties, the dust correction and resulting 
SsFR is heavily affected by the 2 mag Apuy prior (e.g. 
SN 2003ic in Figure [^. Because using this prior to cor¬ 
rect for dust local to the SNIa can have up to a ^1 dex 
effect on the measured Ssfr, we examined the effect of 
changing the Bayesian dust prior to Apuv = 1-0 ± 0.6 
and Apuv = 3.0 ± 0.6. These values span the full range 


of Apjjv in blue galaxies measured by Salim et al. (2007 
see their Figure 13). Changing this prior serves as a way 
to alleviate some of the uncertainty associated with our 
global SFR determination; lowering this prior by 1 mag 
changes ~10 SNe in our sample from laa to lae. 

Following R15, we tried an additional 3 local aperture 
diameters between 2 and 6 kpc because the choice of a 
4 kpc aperture is somewhat arbitrary and other reason¬ 
able choices exist. In part, the FWHM of the FUV PSF 
determines the minimum spatial scale we can probe with 
GALEX, which is approximately 2 kpc at our median 
redshift. However, Figure shows that it is still possi¬ 
ble that a local aperture ■™1 encompass components of 
a galaxy with different star-forming environme nts. The 
highe r -resol ution star formation maps of M33 in Boquien 
et al. (20151 show large Esfr variation on much smaller, 
sub-kpc scales. Nevertheless, we might hope that star- 
formation within a ^few kpc aperture is still much bet¬ 
ter correlated with the SN progenitor environment than 
a global measurement due to the signihcant fraction of 
prom pt progenitors and low velocity dispersions of young 
stars ( de Zeeuw et al.|[l999 ). 

The boundary between SNe laa and lae is also some¬ 
what arbitrary. We used values of lo gfEsFR) between - 


3.1 and -2.7. For direct comparison to Kelly et al. (2015), 
we also examined the boundary between star-iormmg 
and passive of log{T,sFR) = —1.7 and -1.85 (account¬ 
ing for a ~0.4 dex offs et between our SFR measurements 
and Kelly et al. 2015) when discussing Hubble residual 


dispersion. 

Finally, we tried using global rather than local star 
formation (global star formation is a less noisy measure¬ 
ment), and with or without 2.5(T-clipping. Our list of 
analysis variations is given in §4, Table 

4. RESULTS 


We used 179 GALEX-detected SNe from JLA-l-PSl 
and 157 SNe from Rll to measure the LSF step and 
distance dispersion. Although for certain variants of the 
analysis, we see differences between SNe lae and laa at 
the level of ^l-3cr, the evidence for the LSF step is gen¬ 
erally weak. 

Although certain peculiar SNe (e.g. SN 1991bg-like 
and SN 1991T-like) are not explicitly identihed and re¬ 
moved from these samples, the shape and color cuts ap¬ 
plied by JLA and Rll are sufficient to remove many of 
them. However, we make no effort to exclude peculiar 
SNe that JLA/Rll have determined to be cosmologi- 
cally useful so that we can directly assess the affect of 
local SF on the JLA/Rll cosmological analyses. In con¬ 
trast, [Rigaurretari((2M^ and R15 remove identified SN 
19911' explicitly (~3% of their sample). 

In this section, we do not examine the effect of cor¬ 
recting fo r the relationship be tween host mass and SN 
distance (Sullivan et al. 2010| on the LSF step as only 
^15% of our SNe are low-mass hosts (log(M 0 ) < 10; 
R15 similarly found that few H09 SNe are in low-mass 
host s). However, we briefly consider its effect on Hq in 
§5.2[ A complete table with our GALEX measurements 
and Hubble residuals is available online^ with the first 
25 rows given in Table 


4.1. The Local Star Formation Step 

We find a greatly reduced LSF step compared to R15 
for all light curve htters and values of Ry. Using SALT2, 
we hnd an LSF step of 0.000±0.018 mag. With MLGS 
Ry = 2.5 (the value used in the Rll baseline analysis), 
we Hnd 0.029±0.027 mag. However, we do find mild ev¬ 
idence for an offset of 0.059±0.025 mag with Ry = 2.0 
(2.4cr significance). For Ry = 3.1, we found a value of 
0.013±0.030 mag. Our error budget includes systematic 
errors, which we estimated by measuring the standard 
deviation of several variants of our analysis. 

Figured presents our baseline measurement of the LSF 
step and Hubble residual dispersion for SNe la in locally 
passive and locally star-forming environments (SNe lae 
and SNe laa, respectively), with colors indicating the 
probability incorporated in our likelihood model that a 
given SNIa has a locally passive environment, P(Iae). 
We find that 47.2% of Rll SNe in our sample are lae and 
46.0% of JLA-l-PSl SNe in our sample are lae. The over¬ 
all intrinsic dispersion for our full MLGS sample (^0.13- 
0.17 mag; 0.14 for Ry = 2.5) is higher than for SALT2 
( 0.12 mag), likely due to the lack of recent calibration of 
MLGS2k2. Intrinsic dispersion can also be affected by 
the distribution of light curve parameters in the sample 
and the robustness of the photometric measurements. 

We find no significant difference in dispersion between 
SNe laa and SNe lae in SALT2. In the Rll MLGS sam¬ 
ple, however, we find some evidence that SNe lae have 
lower dispersion (cr/ae) than SNe laa. For Ry = 3.1, the 


^ http://www.pha.jhu.edu/~djones/lsfstep.html 
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Figure 4. Our baseline analysis for the JLA+PSl sample (SALT2; upper left), and the Rll sample with different values of Rv (MLCS2k2 
fitter). The color of each SN indicates the probability that it has a locally passive environment, P(Iae). Shaded bars indicate the uncertainty 
on the mean (dark shading; statistical error only) and the standard deviation of the maximum likelihood gaussian (the weighted dispersion; 
light shading). The LSF step is much smaller and has lower significance than the step found by R15, although we detect it at 2.her for the 
Ry = 2.0 case (2.4cr with systematic errors). For Ky = 3.1 and 2.5, we find lower dispersion among SNe in locally passive environments 
than those in locally star-forming environments at 3.3cr and 2.2cr, respectively. For consistency with R15, SNe with only SgFR ^PP^r limits 
are placed at log(EsFR) = -5.3. Systematic uncertainties are estimated from several variants of our analysis (Table [^. 


LSF step is the lowest and ajaa is the highest (0.09 mag 
> a/ae; 2.6a with sys. error). These results disagree 
with R15 at the 3a level. For MLCS with Ry = 2.5, 
(Tjae is ~0.05 mag less than ajaa (1.9a significance). For 
Ry = 2.0 we detected only a ^0.03 mag difference in dis¬ 
persion (1.3a). Our full results for each analysis variant 
are presented in Table 

We found that if we restrict to z > 0.023 (the R15 
minimum z), we see more evidence for the LSF step. Af¬ 
ter this cut, there are 135 SALT2 SNe la and 104 MLCS 
SNe la. The increased significance of these results is ex¬ 
pected because ~3/4 of our MLCS sample is from R15 


when we apply this redshift cut. For MLCS Ry = 2.0, 
2.5, and 3.1 we find LSF steps of 0.086±0.028 (3.1a), 
0.076±0.030 (~50% of R15; 2.5a), and 0.064±0.037 (35% 
of R15; 1.8a). For SALT2, we only find a very small off¬ 
set, 0.017±0.019 (18% of the R15 result) at 0.9a signif¬ 
icance. The MLCS LSF steps are ~50% of those found 
by R15. Except in the case of MLCS with Ry = 2.0, 
the low-z data alone (0.01 < z < 0.023) show slightly 
brighter SNe laa by ~0.02-0.03 mag but with only 0.5a 
significance for MLCS (0.06 mag with 1.4a for SALT2). 
This effect is mostly due to ~5 bright low-z SNe, which 
do not have a large effect on the final result (see the 
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2.5cr-clipping in Table [^. If the peculiar velocity correc¬ 
tions and uncertainties for low-2; SNe were in error, we 
would expect, but do not observe, a significant increase in 
uncertainty-weighted dispersion below z = 0.023 

(we see <0.015 mag difference). We did not find evi¬ 
dence that our highest-2; data (z > 0.07) were having a 
significant effect on our results. 


4.2. Systematic Uncertainties 

Several different variants of our analysis are consis¬ 
tent with the baseline result. The JLA-l-PSl variants 
are shown visually in Figure!^ and the Rll variants are 
shown in Figure]^ For the LSF step, the full results from 
both data sets are presented in Table and our disper¬ 
sion results are presented in Table We have added 
the standard deviation of the LSF step from all variants 
in quadrature to our measured values (giving each type 
of variant, e.g. aperture size, SFR boundary, etc., equal 
weight). Because using the global SFR is not truly a 
local measurement, we have excluded it from our error 
computation but include it in our list of variants for com¬ 
parison. 

For nearly all samples, our most significant detections 
of the LSF step were at a log(SsFR) boundary of -3.1 and 
a 3 kpc aperture radius. For a log(SgFR) boundary of - 
3.1, with SALT2 and MLCS Ry = 2.5 (the most relevant 
versions for cosmology), we detected steps of 0.023±0.019 
and 0.044±0.029, respectively. These are ^25% of R15 
values and insignificant. 

For MLCS with Ry = 2.5 and 3.1, our most signifi¬ 
cant detections came from the variant with 2.ha-clipping. 
They had values of 0.060±0.026 mag (2.3 (t) for Rv = 2.5 
and 0.046±0.028 (l.hcr) for Rv = 3.1. This may mean 
that outliers are affecting our measurement. However, 
we also expect that they affect the Rll Hg measurement 
in the same way, and note that Ry = 2.0 2.5(T-clipping 
has no significant effect. 

The variant with the smallest LSF step was the one 
based only upon global SFR instead of local. However, 
the significance of the difference is only <lcr except in 
the case of Ry = 2.0. The difference may stem from the 
fact that 25% of SNe with globally star-forming environ¬ 
ments in our samples had locally passive environments 
(P(Iae) > 50%). Only 5% of SNe with globally passive 
environments had a >50% probability of bein g locally 
SF. Q u alitat ively, this agrees with Ha data from |Rigault| 
et al. (2013, their Figure 5), who found that globally 
star-forming hosts often had locally passive regions. 

Even after adding the systematic error in quadrature, 
the MLCS Ry = 2.0 LSF step is detected at 2.4cr 
(0.059±0.025 mag). Future cosmology analyses using 
MLCS with low Ry should measure the LSF step in their 
samples to evaluate its effect on cosmology. 

The difference in the dispersion between the two SN 
populations in MLCS is greatest in those same analy¬ 
sis variants discussed above, but as with our baseline 
analysis, we see the opposite effect that R15 found. We 
don’t detect any difference in dispersion for SALT2 with 
the exception of using global instead of local SFR, for 
which we find a 0.05±0.018 mag {2.8a) reduction in dis¬ 
persion for passive hosts. For MLCS Ry = 2.5 and 3.1, 
we find a reduction in dispersion for locally passive SNe 
of ^0.05—0.1 mag (^l-3cr) for a log(EsFR) boundary of 


H09 SNe la with SALT2 distance moduli from JLA (63 SNe) 



Baseline SFR boundary Aper. Size Global SFR 2.5(T-clip 


Figure 5. The systematic error of the SALT2 LSF step estimated 
by the effect of different variants of our analysis on the measure¬ 
ment of the LSF step. Red error bars represent the standard de¬ 
viation of all variants of our analysis added in quadrature to the 
uncertainties from each in dividual varian t . The top panel shows 
only H09 SNe included in|Betoule et al.| ||2014|l, and the bottom 
panel shows our full SN la dataset, me step we detect is ~0.05 
mag (1.3cr) with H09 SNe, but shrinks to <0.01 when we add in 
our full SN la sample. The blue dashed lines and shaded regions 
show the R15 LSF step and Icr uncertainty for SALT2. The results 
from different variants of our analysis are very consistent; our mea¬ 
sured systematic errors are only a small fraction of our statistical 
errors. The global SFR variant is excluded from the systematic 
error calculation, as this is not a local measurement. 

-3.1 and a 3 kpc aperture radius. 

4.3. Consistency with R15 

R15 measured a much larger LSF step of 0.094±0.037 
with SALT2, 0.155±0.041 with MLCS2k2 Ry = 2.5 and 
0.171±0.040 with MLCS2k2 Ry = 3.1. We did not di¬ 
rectly compare to their Ry = 1.7 data, but our Ry = 2.0 
offset is 50% smaller than theirs. Our measured SALT2 
LSF step has a 2.3cr discrepancy with the R15 measure¬ 
ment, our MLCS2k2 Ry = 2.5 LSF step has a 2.6ct dis¬ 
crepancy, and our MLCS2k2 Ry = 3.1 LSF step has a 
3.2(7 discrepancy. 

Table demonstrates the step by step impact of 
changes in the R15 analysis or data, showing the effects 
of using the JLA-l-PSl and Rll light curve cuts, the 
JLA-l-PSl and Rll distance moduli (with an updated 
SALT2 light curve fitter for JLA-l-PSl), our improved 
log(EsFR) measurements, and using a larger SN la sam¬ 
ple (with and without the R15 z > 0.023 cut). 

Updated distance moduli greatly decrease the signif¬ 
icance of the LSF step in JLA-l-PSl data in SALT2, a 
50% reduction (a change in significance of 2.4 (t). The ver¬ 
sion of SALT2 used in recent analyses has an improved 
SN la model and uncertainty propagation, a larger train¬ 
ing sample, and an updated value for /3. Rll distances 
are nearly identical to H09 distances, so using these has 
no significant effect on the LSF step. 

Using our Esfr measurements increases the signifi¬ 
cance of the LSF step by O.Scr for SALT2 and reduces 
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Table 3 

The Effect of Step-by-Step Changes in R15 Data, Distances, SFR Measurements, and Sample Cuts 


Measurements SALT2 MLCS Rv=2.5 

SN SN 


sample 

l-^resid 

SsFR 

cuts 

SNe 


Sig. 


Sig. 

SNe 


Sig. 

'^laa 

Sig. 

H09 

H09 

R15 

H09 

77 

0.093±0.026 

3.5o- 

-0.034±0.073 

-0.5(7 

81 

0.169±0.026 

6.5(7 

0.057±0.033 

1.7(7 

H09 

H09 

R15 

JPR^H09 

59 

0.129±0.030 

4.3o- 

0.012±0.047 

0.2(7 

74 

0.144±0.025 

5.61T 

0.038±0.034 

1.1(7 

H09 

JPR 

R15 

JPR,H09 

59 

0.062±0.032 

1.9o- 

0.030±0.031 

1.0(7 

74 

0.149±0.025 

5 . 9 (t 

0.023±0.031 

0.7(7 

H09 

JPR 

Here‘S 

JPR,H09 

59 

0.071±0.033 

2.2o- 

0.009±0.031 

0.3(7 

74 

0.119±0.026 

4.5(7 

-0.010±0.030 

-0.3(7 

H09 

JPR 

Here 

JPR 

63 

0.045±0.033 

1.3o- 

0.015±0.030 

0.5(j 

78 

0.097±0.027 

3.61T 

-0.029±0.030 

-1.0(7 

JPR® 

JPR 

Here 

JPR,2 > 

0.023 

135 

0.017±0.019 

0.9o- 

-0.020±0.019 

-1.1(7 

103 

0.076±0.029 

2.6(7 

-0.041±0.029 

-1.4cr 

JPR® 

JPR 

Here 

JPR 

179 

O.OOOiO.018 

O.Oo- 

-0.013±0.018 

-0.7(7 

156 

0.029±0.025 

1.2(7 

-0.053±0.024 

-2.2(7 


Note. — We show the difference between our analysis and R15 by improving one element of the analysis at a time. We start with the R15 results and 
sequentially show the effect of adding light curve cuts from JLA+PSl/Rll, using JLA/Rll distance moduli, using our updated SFR measurements, using 
only JLA/Rll (not H09) light curve cuts, and finally adding in the full SN samples with and without the R15 redshift cut of 2 > 0.023. The biggest 
differences come from adding the full sample for both SALT2 and MLCS and using improved SALT2 distance moduli. The Rll SN light 
curve cuts also make a Icr difference in the MLCS results. For consistency, we have used the likelihood minimizer used in the rest of this study to reproduce 
the R15 results (The SciPy Optimize package). This minimizer returns smaller uncertainties than Minuit, which was used in R15, but we find negligible 
differences in the maximum likelihood values themselves. The difference in LSF step we find for R15 data with MLCS (our value is 0.014 mag higher) is 
because we adopt two separate dispersions for SNelaa and SNelae whereas R15 use a single value for the full sample. 

^ denotes the magnitude of the LSF step. 

^ The differ ence in uncertainty-weighted dispersion between SNelae and laa (using the standard deviation of the maximum likelihood gaussians; CTg and (Ta 
in Equation [AR. 

for SALT2, Rll for Ry ^ 2.5. 

^ Measurements of Ssfr from this work (see §3). 

® The full JLA+PSl (SALT2) and Rll (MLCS) SN samples. 



Figure 6. The systematic error of the MLCS LSF step estimated 
by different variants of our analysis for Ry = 2.0, 2.5, and 3.1 
in the Rll SNIa sample. The LSF step has 2Aa significance for 
Ry = 2.0. The baseline analysis used to determine Hq uses Ry = 
2.5, for which we see a small LSF step at l.lo- significance. We 
see <lc7 significance for Rv^ = 3.1. The blue dashed lines and 
shaded regions show the R15 LSF step and la uncertainties for 
MLCS2k2. The global SFR variant is excluded from the systematic 
error calculation, as this is not a local measurement. 


it by 1.4cr (~20%) for MLCS. Between our data and the 
R15 data, there is significant scatter in probability for 
10% < P(Iae) < 90%, in large part due to our modest 
changes in dust correction methodology. However, we 
find only 3% me dian offset in P(Iae) between our data 
and R15 and in §4.6| we find that our method of Ssfr 
measurement has little impact on the final results. Our 
full set of Esfr measurements can be compared to R15 
using the data we provide online and in Table 

There are 4 SNe in Rll and 4 SNe in JLA ^at pass 
Rll/JLA light curve cuts but do not pass H09 cuts (SNe 
1992j, 1993h, 1999aw, 2001ic, 2006bd, 2006gt, 2007ba, 
and 2007cg). We found that including them reduces the 
SALT2 LSF step by a significant 37% (0.9 ct) and reduces 
the MLCS LSF step by ^15% (0.9cr). When applying 
any LSF-dependent effect to cosmology, it is appropriate 
to match the cuts used in the cosmological analysis to 
those used in the measurement. 

For both the LSF step and the dispersion in MLCS, 
there is a > Icr change when we use the full SN la sample. 
Although the total statistical change from 3.6cr to 1.2cr 
is large, we do not expect this to be a result of peculiar 
velocity bias from our low-z data. Some of the change 
may result from a greater sample dispersion, which re¬ 
duces the significance of small offsets. A dispersion term 
is typically added in quadrature to distance modulus un- 
certainties in cosmol ogical analyses, including Rll and 


Betoule et al. 

02014 

1, and has the same effect. In ad- 

dition, i'able 

does not incorporate systematic error. 


which may have an impact; high- 2 ; data effectively have 
a larger aperture size due to a PSF width that is greater 
fraction of the 4 kpc aperture diameter. Figure Iw shows 
that aperture variations may have up to a Icr 5iect on 
the measured LSF step, and to expand our sample size 
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we have preferentially added low-z data with smaller ef¬ 
fective apertures (0.01 < z < 0.023). 

Table 1^ shows that the MLCS increase in laa disper¬ 
sion is mostly caused by the addition of new SNe rather 
than to our Ssfr measurements or new distance mod¬ 
uli. The surveys that comprise our sample typically have 
larger dispersion than H09, which reduces the signih- 
cance of the H09 sample. There are a number of possible 
sources for increased dispersion of a SN la sample, includ¬ 
ing underestimating photometric difference image uncer¬ 
tainties near bright hosts and nig htly or absolute pho - 
tometric calibration uncertainties (Scolnic et al.||2014a). 
For MLCS, Rll may also have higner sample dispersion 
because they make no cut on the of the MLCS light 
curve fits, while H09 remove SNe with reduced > 1.5. 


4.4. The Effect of MLCS Sample Cuts 

In MLCS, the total difference of ~0.14 mag between 
our analysis and R15 may appear surprising, but in ad¬ 
dition to the possible reasons discussed above, much of 
the change between the R15 measurement and ours ap¬ 
pears to arise from the different demographics of the two 
samples and the peculiarities of the MLCS light curve 
fitter. H09 find that for both high-Ay SNe and high- 
A SNe, MLCS tends to overcorrect leading to negative 
residuals, and these negative residuals are not subtle. In 
our Ry = 2.5 sample, SNe with Ay > 0.5 have a mean 
residual of -0.22 mag, which has been seen elsewhere as 
evidence for a lower Ry in high extinction environments. 
Likewise, SNe with A > 0.7, where the relation between 
light curve shape and luminosity becomes non-linear and 
is poorly sampled especially when MLCS2k2 was trained, 
have a mean residual of -0.23 mag. Accordingly, the bal¬ 
ance of rare high Ay SNe to rare high A SNe can affect 
an apparent LSF step as the frequency of these objects 
correlates with host properties. 

Passive hosts have preferentially higher A than SF 
hosts (H09, their Figure 19), while SF hosts have pref¬ 
erentially higher Ay. In R15, the H09 data that have 
GALEX imaging and pass their cuts contain several SNe 
with large A but only two SNe with Ay > 0.45 for 
Ry = 1.7 (for Ry = 3.1, only two SNe with Ay > 0.7). 
Therefore a sample like R15 without high-Ay hosts but 
with high-A hosts will have brighter passive SNe la on 
average, producing a larger apparent LSF step. 

One approach to decrease sensitivity to MLCS Hubble 
residual trends is to first remove the trends, and then 
determine the LSF step. In Figure we ht a simple 
linear model to MLCS Hubble residuals as a function 
of A and Ay, using Rll SNe in H09 (with Ay < 1.5 
and A < 1.5 to match H09). When we correct for those 
slopes, we see that the measured SF step using Rll SNe 
in H09 shrinks by a factor of 2.5 and is reduced from 4.5cr 
to 1.6cr significance. 

SALT2 does not have the strong residual trends with 
Xi and C that MLCS does with Ay and A, and we 
also find that restricting our sample to the H09 “best” 
SALT2 cuts (—0.1 < C < 0.2) do es n ot introduce an 
LSF step (but changing (3 may; see (5.1). However, it is 
likely that recent substantial improvements to the SALT2 
model have removed some of the biases in its derived 
distances. Due to the lower dispersion of SALT2-fit SNe, 
the lack of these residual trends, and because MLCS fits 
assume an extinction law, it is likely that SALT2 is more 



A Ay 


R15 SNe in Rll without modei correction; MLCS R^. = 2.5 



logics™ h'l. vr ' kic b) 


R15 SNe in Rll with model correction; MLCS R,, = 2.5 



[1/q yr * Lpc ■ ]) 


Figure 7. A simple linear correction for Hubble residual trends 
in MLCS reduces the significance of the R15 LSF step. In the top 
panels, we show MLCS A and Ay fit to Rll SNe. In the middle 
panel, we show our measured SF bias using Rll SNe in H09. In the 
bottom panel, we make a linear correction for the MLCS Hubble 
residual trends, and the LSF step is reduced from 4.5fT to l.ficr 
significance. Colors indicate P(Iae), with P(Iae)~100% in red and 
P(lae)~0 in blue. 


effective at standardizing SNe la. 

In a future update of MLCS using a larger training 
sample, it would be important to verify that these trends 


with host, Ay and A are diminished. 

4.5. Kelly et al. 

(2015) Scatter 

Using MLCS, Kelly et al. 

(2015) see reduced Hubble 

residual scatter of only 3.57 

in distance in highly star- 


forming regions (log(EsFR) > -2.1 and log(SsFR) > - 
2.25). Due to differences in methodology, there is a ~0.4 
dex offset in measurements between our data and 


Kelly et al. (2015). Because of this, we adopt log(EsFR) 
> -1.7 and log(EsFR) > -1.85 as our T,sfr boundaries 
for comparison. 


In part, the low scatter seen by Kelly et al. (2015) is 
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because they explicitly remove SNe with Hubble diagram 
residuals >0.3 mag (>15% in distance). Because of this 
and because the Rll sample does not cut SNe with high 
extinction or large A, our unweighted standard deviation 
is a significantly larger ^0.25 (12% in distance) for the 
Rll sample at log(I]sFR) > -1-7 and log(EsFR) > -1.85. 
For SALT2, the standard deviation is a slightly lower 
0.20 mag, or 10% in distance, with no difference between 
SNe in locally passive/locally star-forming environments. 

We also see no difference in uncertainty-weighted dis¬ 
persion for these SgpR boundaries in SALT2, and we 
find that the dispersion for SNe in both passive and star¬ 
forming environments in SALT2 data is smaller than the 
lowest dispersions we observe with MLCS. The scatte r 
in our sample is much higher than in Kelly et al. (2015), 
and we find a <0.02 mag (^14%; ~0.l-0.5tT) reduction 
in dispersion for MLCS with Ry = 2.0. SNe in star¬ 
forming environments have higher dispersion with low 
significance for MLCS Rv = 2.5. For Ry = 3.1, SNe in 
star-forming environments have ~0.07 mag higher dis¬ 
persion at ~1 ct significance. A summary of our intrinsic 
dispersion measurements are in Table 

If we apply H09 A and Ay cuts to our dat a, we still 
see the opposite e ffect as [Kelly et a l.| (|2015|). We can 
only reproduce the Kelly et al. 12015) results using their 
strict A and Ay cuts, which have not been used in any 
cosmological analysis to date. However, these cuts may 
prove useful in the future if this low-scatter population 
persists when additional SNe are added to the data. 


4.6. Additional Consistency Checks 

We performed several consistency checks to verify that 
individual SN datasets and differences between our anal¬ 
ysis and R15 did not bias our results. First, we re¬ 
moved SNe discovered prior to the year 2000, leaving 130 
SNe from JLA/PSl and 116 SNe from Rll. Our results 
were consistent with our baseline analysis; we measured 
a SALT2 LSF step of 0.010±0.025 mag and an MLCS 
Ry = 2.5 step of 0.040±0.031 mag. The Ry = 2.0 step 
was a slightly higher, but consistent, 0.079±0.030 mag 
(2.7(t). The dispersion of SNe in highly SF regions was 
not significantly reduced. 

Second, the photometry and calibration from low-z 
surveys is not as robust as recent data from SDSS and 
PSl. The JLA/PSl sample has 37 SNe with redshifts less 
than 0.1 that have GALEX data and pass our cuts, while 
the Rll sample includes no SDSS/PSl SNe as it predates 
them. For comparison, we fit SDSS and PSl SNe with 
MLCS to see if the LSF step derived from these surveys 
alone are consistent with the Rll results. With SALT2, 
we find an LSF step of 0.034±0.028 mag with lower SF 
dispersion by 0.049±0.024 mag (2.0(t). With MLCS, we 
find a large LSF step with 35 SNe of ~0.14±0.055 mag 
with 1.7-2.9cr significance. As the sample consists of only 
~10-15 locally passive SNe, this step could still be caused 
by low statistics or a limited range of light curve param¬ 
eters comprising the sample. As discussed in §4.4[ the 
trends MLCS residuals have with different light curve 
parameters may be a factor, as the size and significance 
of the LSF step is somewhat reduced when this sample 
is restricted to low A and Ay. This step is also unlikely 
to affect recent cosmological analyses, which are based 
on SALT2 or comprised mainly of low-z data (e.g. Rll, 
H09). However, it is an interesting result that should be 


explored further with photometric PSl SNe and future 
DES data. This sample is too small at log(S sFR ) > -1.85 
for a reliable check on our Kelly et al. (20151 comparison. 

If we make a host galaxy inclination cut at >80° follow¬ 
ing R15 (instead of our more conservative cut of >70°), 
the results are consistent with our baseline result, with 
MLCS LSF steps ranging from 0.00 mag (Ry = 3.1) to 
0.045 mag (Ry = 2.0) with uncertainties ~0.025 mag. 
The SALT2 LSF step is -0.016 mag (<la significance). 

Finally, we apply a dust correction to the FUV flux 
from all SN regions in star-forming hosts when deter¬ 
mining EgpR, now including t he 20 Rll SNe and 25 
JLA/PSl SNe with R > 3 (see (j3.1.2|). We again find a 
comparable result; the SALT2 LSF step is 0.012±0.019 
mag, and the MLCS Ry = 2.5 LSF step is 0.040±0.028 
mag. 


5. DISCUSSION 

We find that local star formation has at little to no 
effect on SNIa distances in the Rll and JLA-l-PSl sam¬ 
ples. Our results have several important implications for 
cosmological analyses, Hq, and future measurements of 
relationships between SNe la and their host galaxy prop¬ 
erties. 


5.1. The Effect of fi and Ry on SNIa Distances 

Although the modest differences we observe in mean 
magnitude and dispersion for MLCS with certain values 
of Ry could be due to the relation between SN la pro¬ 
genitor properties and derived distances, we consider it 
much more likely that host galaxy extinction, which is 
highly correlated with star formation, is causing any ob- 
serve d bias. We propo se that some of the effects seen in 
R15, Kelly et al. (2015 ), and our data may be due to dust 
rather th an to a secondary eff ect such as the progenitor 
age (e.g. Childress et al. 2014). 

With MLCS, the LSI step we found is 0.046±0.039 
mag higher assuming Ry = 2.0 than assuming Ry = 3.1 
(systematic errors added). The Ry = 2.0 LSF dispersion 
is 0.053±0.044 (stat-|-sys) mag lower than Ry = 3.1. It 
has b een observed by several groups (e.g. [Burns et al.j 
2014) that SNe la in high-extinction environments have 
lower values of Ry. Because of this, it seems likely that 
the Ry = 3.1 extinction law is failing to properly correct 
for the dust in some star-forming regions. 

For SALT2, our value of /3 has a value ~0.6 higher in 
the latest cosmological analyses than the value found in 
H09. This can have an important effect on the measured 
LSF step. For example, a SNe la in a locally star-forming 
environment with ~0.17 magnitudes of Ay, would have 
its corrected magnitude shifted by 0.1 mag with this new 
value of /3. For comparison, R15 SNe with locally star¬ 
forming environments have a mean fitted Ay = 0.25 for 
Ry = 3.1 and Ay = 0.22 for Ry = 1.7. We don’t see 
such a large effect in our data, and would not expect 
/3 to have the exact effect of Ry, but we do find that 
using a lowered /3 of 2.5 (the value used in H09) in our 
analysis raises the SALT2 LSF step to 0.024±0.018 (1.3(7 
significance). 

In future cosmological analyses, it may be possible to 
separate star-forming and passive hosts and fit for two 
different values of /3 or Ry. This could reduce scatter 
and provide more precise SNIa distances for subsets of 
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the population, provided the systematic uncertainties in 
such an analysis are well-understood. 

The SALT2 light curve fitter shows the least differ¬ 
ence between SNelae and SNelaa and also has 

the lowest dispersion in both star-forming and passive 
regions. The lowest dispersion we find using MLCS is 
still higher than the SALT2 dispersion for both SNe lae 
and laa. For this reason, SALT2 may be a more reliable 
light curve fitter for cosmological analyses. In its current 
version, MLCS fails to standardize SNe la to the extent 
that SALT2 does and has fitter biases that correlate with 
host properties (such as Hubble residual nonlinearities 
with high A and an assumed value for Ry). Perhaps a 
re-trained version of MLCS t hat incorporates term s such 
as random SN color scatter ( Scolnic et al.||2014b ) would 
reduce the MLCS outlier fraction and provide more pre¬ 
cise distances. 

5.2. The Effect on Measuring Hq 

Because our final measurement of the LSF step with 
Ry = 2.5 is only a 1.2cr detection, there are no grounds 
in the Bayesian sense to correct Hq for the LSF step. 
However, a useful test of systematic uncertainties in the 
future will be to use only star-forming hosts in the Hub¬ 
ble flow sample, which have similar physical properties 
to the nearby Cepheid-calibrated sample and will better 
control for unknown biases in metallicity, dust, or pro¬ 
genitor age. 

Adopting the 47.2% SN lae fraction we find for Rll and 
the 7.0% SNIae fraction found by R15 for the Cepheid 
sample with Equation we find no evidence for a re¬ 
duced value of Hq. Following R15, if we were to replace 
the host mass step with the LSF step, our measurement 
suggests a 0.1% increase in Hg because the size of the 
LSF correction is slightly less than the size of the host 
mass correction. 

One caveat is that Rll added the MLCS intrinsic SN la 
dispersion but not the full apparent intrinsic dispersion 
in quadrature to the distance modulus uncertainties in 
their Hubble flow SNe. We find that forcing our max¬ 
imum likelihood gaussian model to use only the MLCS 
intrinsic dispersion of 0.08 mag raises the magnitude of 
the Ry = 2.5 LSF step we derive to 0.045±0.019 (a 2.4cr 
detection, but 2.1cr with systematic uncertainty added). 
This could be because it allows outliers to have a greater 
effect on the measurement. However, applying this cor¬ 
rection after removing the host mass step still only re¬ 
sults in a reduction in Hg of 0.11 km s“^ Mpc“^. The Rll 
value for Hg is within the Icr uncertainty of the LSF step. 
The highest LSF step we are able to find using all our 
analysis variants with 0.08 mag dispersion is 0.066±0.22 
mag (the 2.5CT-clipped variant), and even this extreme 
measurement lowers Hg by only 0.4 km s“^ Mpc“^. 

Finally, if we measure the LSF step after ho st mass 
correction using masses from Neill et al. (2009 53% of 
the Rll sample) and again using a dispersion of 0.08 
mag, we find a LSF step of 0.023±0.027 (stat-|-sys) mag 
for Ry = 2.5. This results in a small reduction of 0.3 km 
s“^ Mpc“^. Because we detect this effect at <lcr (with 
systematic error added in quadrature), we do not believe 
a correction is justified. 

5.3. Future Measurement of the LSF Step 


Although we have only detected the LSF step at low 
significance with GALEX FUV data, GALEX alone is 
not the best tool for studying local regions due to its large 
PSF width and the uncertain UV extinction correction. 
The LSF step would be best identified in local Ha (e.g. 
Rigault et al. 2013), high-resolution UV data from the 
HubbLe Space Telescope (HST), or local SED fitting. 

Table shows that sample selection has a significant 
effect on our results. We suggest that studies examin¬ 
ing host galaxy effects use the same SNIa samples and 
selection criteria as the latest cosmology analyses when 
possible. It may be possible to detect the LSF step or dif¬ 
ferences in dispersion at higher significance using differ¬ 
ent light curve or distance modulus cuts, but the results 
of such analyses would not necessarily apply to typical 
measurements of cosmological parameters. 

Local SED fitting may be the optimal approach for 
studying the relation between host galaxy properties and 
SNIa distances, as it can put simultaneous (albeit some¬ 
times degenerate) constraints on a number of parameters 
that may correlate with SNIa distances such as stellar 
age, extinction, star formation history, and mass con¬ 
tained in a local region. Approaches that don’t depend 
entirely on GALEX data will also be able to measure lo¬ 
cal regions at higher redshifts and put better contraints 
on possible redshift-dependent biases. 

The size of the samples with which we can examine 
the effects of host galaxy properties on SNIa corrected 
magnitudes will increase dramatically in the next few 
years. The PSl photometric sample alone will consist 
of up to ~2,000 SNe la with cosmologically-useful light 
curves. The Dark Energy Survey (DES) will contribute 
thousands more up to redshifts of ~1. Although mea¬ 
surements of local regions become more difficult at high- 
z, a ground-based optical survey with PSF FWHM ~1 
arcsec will be able to use a much larger SN sample pro¬ 
vided the absence of UV data is not prohibitive. Surveys 
such as PSl or DES are able to examine local regions of 
5 kpc diameter, similar in size to the apertures used in 
this study, up to z ~ 0.35. 

6. CONCLUSIONS 

Analyzing the same SNe la used to determine the most 
recent values of w and Hg, we find little evidence for a 
LSF step, which suggests that correcting cosmological 
parameters for this effect is not necessary. There is only 
l.lcr evidence for the LSF step in Rll MLCS data assum¬ 
ing Ry = 2.5 (the Ry Rll used in their baseline analysis) 
and O.Ocr evidence for the LSF step in JLA-l-PSl SALT2 
data. Our most significant detection uses MLCS data as¬ 
suming Ry = 2.0, for which we find 2Aa evidence for a 
step. The sizes of both of these steps are greatly reduced 
compared to the measurement of R15. Lower values of 
P in SALT2 and Ry in MLCS may increase the size and 
the significance of the LSF step. 

Compared to R15, differences in our Tisfr measure¬ 
ment and dust correction technique reduced the size of 
the MLCS LSF step by ~20% and increased the SALT2 
LSF step by ~15%. Using MLCS sample cuts from Rll 
reduced the offset by an additional ~20% and adding the 
full Rll sample reduced the offset to 0.029±0.027 mag, 
likely due to the higher dispersion and better statistics 
of the full sample. Using new distance m oduli and sam¬ 
ple cuts from only Betoule et al. (2014) (and not H09) 
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reduced the SALT2 LSF step by 60% and using the full 
JLA+PSl sample reduced the SALT2 step to a value of 
0.000±0.018 mag. 

MLCS sample cuts have a significant impact on the 
results. MLCS Hubble diagram residuals are more nega¬ 
tive at greater Ay and A, which must be carefully taken 
into account in cosmological analyses. In particular, pas¬ 
sive hosts are known to have preferentially higher A but 
lower Ay (H09). We suspect that because the R15 sam¬ 
ple had few high-Ay SNe but a wide range of A, their 
locally star-forming SNe had preferentially fainter Hub¬ 
ble residuals. 

We found that JLA-fPSl SNe fit with SALT2 had 
lower dispersion than MLCS-fit Rll SNe in star-forming 
or passive environments. We also found that locally star¬ 
forming SNe in our sample did not have lower dispersion 
at log(EsFR) > -2.9. In MLCS with Ry = 3.1, SNe la in 
locally passive environments have lower dispersion than 
those in locally star-forming environments by ^0.09 mag, 
a 2.5o- result. Using MLCS with Ry = 2.5, we see a 
0.053±0.029 mag difference. 


The lowest SNIa dispersions come from using S ALT2 


distance moduli. In contrast to 

Kelly et al. 

(2015 

), with 

MLCJS we found no evidence that yiNe in 

highly star- 


forming environments have lower dispersion than locally 
passive SNe using Ry = 2.0. With Ry = 3.1 we found 
that SNe in star-forming environments had greater dis¬ 


persion (' ^1-2(7 significance) , but note that we did not 
make the Kell y et al. ( 201 5 ) sam ple cuts. We can only 
reproduce the Kelly et al. 12015 1 results by using their 
strict cuts on the SIN light curve parameters A and 
Ay and removing SNe with Hubble residuals >0.3 mag, 
which restr i cts ou r sample to largely the same data as 


Kelly et al. (2015). 

I'he LSI' step may also be difficult to detect because of 
the large PSF width of GALEX and it may also be that 
the LSF step is only apparent in analyses with certain 
types of light curve selection or outlier rejection. Future 
studies with local Ha, SED fitting, or HST UV observa¬ 
tions will have an improved ability to detect local effects. 
Our results also show that certain SN sample cuts may 
inadvertently increase biases in cosmology. We expect 
that with the large SNIa samples from PS I and DES 
that will be published in the next few years, the system¬ 
atic uncertainties on Hg and the dark energy equation of 
state will come into clearer focus. 


This work would not have been possible without com¬ 
ments, suggestions, and other assistance from Mickael 
Rigault. We would also like to thank the anonymous ref¬ 
eree and Pat Kelly for many useful comments and sug¬ 
gestions. 


APPENDIX 

CALCULATION OF PROBABILITIES AND MAXIMUM LIKELIHOOD ESTIMATION 


The only significant difference between our method of measuring the maximum likelihood LSF step and Hubble 
residual dispersions and the RI5 method is that we allowed the intrinsic dispersion of both SNIa populations (lae and 
laa) to be ht by our maximum likelihood model. We describe our full procedure below. 

We first converted the dust-corrected FUV flux into SgpR following RI5 (their Equation I). We set the boundary 
between the locally star-forming and locally passive population at log(E 5 Fi^) = —2.9 as in RI5, and measured the 
probability that the SNIa exploded in a locally passive environment based on the full probability distribution from 
our dust-corrected photometric measurements. 

We used these probabilities to construct a maximum likelihood model assuming two gaussian populations of SNe 
with different mean Hubble residuals and dispersions. The likelihood is determined by the equation: 


Ci = P{Iaa) X 


+P{Ia€) X 


I 

v/27r(a2 + a2)®^P^ 2(af + 

I . 

V27r(af + a2)"^P^ 2(af+a?) ^ 


(Al) 


where M’j^Y is the corrected magnitude and ai is the corrected magnitude uncertainty of a given SNIa. P(Iaa) and 
P(Iae) are the probabilities that the SN environment is locally star-forming or locally passive, respectively. /Tq,, 
and (Te are free parameters equal to the means and standard deviations of the normal distributions of SNe laa and lae. 
To determine what these parameters are, we found the maximum likelihood model by minimizing: 


N 

iog(/:) = -2 ^ iog(/:i) (A2) 

i=l 

where N is the number of SNe la in the sample. 

Instead of adding an intrinsic dispersion term in quadrature to the Hubble residuals such that the reduced of 
the sample is I, as is commonly done in cosmological analyses (and in RI5), we fit to the standard deviations of our 
gaussian maximum likelihood model for SNe laa and lae. We verified that allowing the dispersion to be fit by our 
model instead of specifying it beforehand does not affect our results. 


REFERENCES 

Aldering, G., Adam, G., Antilogus, P., et al. 2002, in Society of Photo-Optical Instrumentation Engineers (SPIE) Gonference Series, Vol. 
4836, Survey and Other Telescope Technologies and Discoveries, ed. J. A. Tyson S. Wolff, 61-72 




























15 


Bertin, E., Arnouts, S. 1996, A&AS, 117, 393 
Betoule, M., Kessler, R., Guy, J., et al. 2014, A&A, 568, A22 
Boquien, M., Calzetti, D., Aalto, S., et al. 2015, ArXiv e-prints 
Burns, C. R., Stritzinger, M., Phillips, M. M., et al. 2014, ApJ, 789, 32 
Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245 
Childress, M., Aldering, G., Antilogus, P., et al. 2013, ApJ, 770, 108 
Childress, M. J., Wolf, C., & Zahid, H. J. 2014, MNRAS, 445, 1898 
Conley, A., Guy, J., Sullivan, M., et al. 2011, ApJS, 192, 1 

de Zeeuw, P. T., Hoogerwerf, R., de Bruijne, J. H. J., Brown, A. G. A., &: Blaauw, A. 1999, AJ, 117, 354 

Dominguez, I., Hoflich, P., &: Straniero, O. 2001, ApJ, 557, 279 

Ganeshalingam, M., Li, W., Filippenko, A. V., et al. 2010, ApJS, 190, 418 

Guy, J., Astier, P., Baumont, S., et al. 2007, A&A, 466, 11 

Guy, J., Sullivan, M., Conley, A., et al. 2010, A&A, 523, A7 

Hamuy, M., Phillips, M. M., Suntzeff, N. B., et al. 1996, AJ, 112, 2398 

Hamuy, M., Folatelli, G., Morrell, N. L, et al. 2006, PASP, 118, 2 

Hayden, B. T., Gupta, R. R., Garnavich, P. M., et al. 2013, ApJ, 764, 191 

Hicken, M., Wood-Vasey, W. M., Blondin, S., et al. 2009a, ApJ, 700, 1097 

Hicken, M., Challis, P., Jha, S., et al. 2009b, ApJ, 700, 331 

Hicken, M., Challis, P., Kirshner, R. P., et al. 2012, ApJS, 200, 12 

Hsiao, E. Y., Conley, A., Howell, D. A., et al. 2007, ApJ, 663, 1187 

Hudson, M. J., Smith, R. J., Lucey, J. R., & Branchini, E. 2004, MNRAS, 352, 61 

Jha, S., Riess, A. G., & Kirshner, R. P. 2007, ApJ, 659, 122 

Jha, S., Kirshner, R. P., Challis, P., et al. 2006, AJ, 131, 527 

Johansson, J., Thomas, D., Pforr, J., et al. 2013, MNRAS, 435, 1680 

Kelly, P. L., Filippenko, A. V., Burke, D. L., et al. 2015, Science, 347, 1459 

Kelly, P. L., Hicken, M., Burke, D. L., Mandel, K. S., & Kirshner, R. P. 2010, ApJ, 715, 743 

Kessler, R., Becker, A. C., Cinabro, D., et al. 2009a, ApJS, 185, 32 

Kessler, R., Bernstein, J. P., Cinabro, D., et al. 2009b, PASP, 121, 1028 

Lampeitl, H., Smith, M., Nichol, R. C., et al. 2010, ApJ, 722, 566 

Li, W., Leaman, J., Chornock, R., et al. 2011, MNRAS, 412, 1441 

Maoz, D., Mannucci, F., Nelemans, G. 2014, ARA&A, 52, 107 

Neill, J. D., Hudson, M. J., & Conley, A. 2007, ApJ, 661, L123 

Neill, J. D., Sullivan, M., Howell, D. A., et al. 2009, ApJ, 707, 1449 

Pike, R. W., & Hudson, M. J. 2005, ApJ, 635, 11 

Rest, A., Scolnic, D., Foley, R. J., et al. 2014, ApJ, 795, 44 

Riess, A. G., Press, W. H., & Kirshner, R. P. 1996, ApJ, 473, 88 

Riess, A. G., Kirshner, R. P., Schmidt, B. P., et al. 1999, AJ, 117, 707 

Riess, A. G., Macri, L., Casertano, S., et al. 2011, ApJ, 730, 119 

Rigault, M., Copin, Y., Aldering, G., et al. 2013, A&A, 560, A66 

Rigault, M., Aldering, G., Kowalski, M., et al. 2015, ApJ, 802, 20 

Salim, S., Rich, R. M., Chariot, S., et al. 2007, ApJS, 173, 267 

Schlafly, E. F., k Finkbeiner, D. P. 2011, ApJ, 737, 103 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 

Scolnic, D., Rest, A., Riess, A., et al. 2014a, ApJ, 795, 45 

Scolnic, D. M., Riess, A. G., Foley, R. J., et al. 2014b, ApJ, 780, 37 

Stritzinger, M. D., Phillips, M. M., Boldt, L. N., et al. 2011, AJ, 142, 156 

Sullivan, M., Le Borgne, D., Pritchet, C. J., et al. 2006, ApJ, 648, 868 

Sullivan, M., Conley, A., Howell, D. A., et al. 2010, MNRAS, 406, 782 

Tripp, R. 1998, A&A, 331, 815 

Tully, R. B., & Fisher, J. R. 1977, A&;A, 54, 661 



16 


Bh 


a 


^ ft 

m CD 
^ CO 

rfS fe 

n CO 


CD 



























































fl 










0 












0 


o 



























cS 

oT 



98 

o 

100 

100 

o 

97 

15 

100 


o 

12 

100 

- 

CM 

lO 

96 

100 

05 

85 

O 

O 

11 

CM 

30 




I-H CD 

M* O 

C55 O 

■':t< CD 

CQ « 

^ CQ 

^ CQ 


O ID 

CQ ID 

00 05 




CQ CD 

t- ^ 


O CQ 

I-H CD 


CD (CQ 

ID CO 

'If d 

CO CO 




GO 

ID ID 

ID CO 

05 CQ 

CO ID 

ID CQ 

00 05 


t- 05 

-H O 

O CQ 


ID ID 

h- IS 

05 O 

00 CD 


o CQ 

O CO 

CQ CQ 

00 1- 

CQ 



Csi 



o o 

O O 

(N •If 

.-H CQ 

o o 

CQ ■'t 

CD ID 


.-H CQ 

T-l I—I 

CQ CD 


I-H CQ 

i-H CQ 

CQ •If 

.-H CD 


CD 'If 

CQ CD 

o o 

o o 

CD ID 

CQ CD 

CQ CD 



6 6 

6 6 

d d 

d d 

d d 

d d 

d d 

zo 

o o 

d d 

O O 

05 

o o 

d d 

d d 

o o 

CM 

o o 

6 6 

6 6 

6 6 

6 6 

CD (d 

(d o 




+ 1 

+ 1 

+ 1 

+ 1 

+ 1 

+ 1 

+ 1 

+ I 

+ I 

+ 1 

+ 1 

+ 1 

+ 1 

+ I 

+ 1 

+ 1 

+ I 

+ 1 

+ 1 

+ 1 

+ 1 

w 


o 

CM 


00 

CM 



CO 


CO 

b- 



o 

X 

X 

o 


ZO 

X 

X 

X 

X 

X 

ZO 


t- 

CO 

05 


o 

CM 

X 


CO 


lO 


o 

i> 

X 

CM 


lO 

t- 

CM 

X 

b- 

05 

CO 

hD 



O 

CM 

00 

r- 

lO 

CO 

(35 

V 

1^ 

X 


V 

(35 

(35 

X 

CM 

V 

1—1 

o 

lO 

I-H 

o 

m 

b- 

q 


© 

CO 

CM 

ft 

cd 

ft 

CO 

ft 

CM 

ft 

CM 

ft 

1—1 

1—1 

cd 

CM 

cd 

CM 

CM 

CM 

ft 

CM ' 




1 

1 

1 

1 

1 

1 

1 


I 

I 

1 


1 

1 

1 

I 


1 

1 

I 

1 

1 

1 

1 

0 



























O 

4J 



X 


X 

x 

X 

X 


X 

ft 

X 




ft 

ft 

2; 

2 


2 

2 




2 

Q 



t- 

00 

o 

CO 

X 





lO 

o 

lO 

X 

(35 

b- 

X 

t- 

lO 

CM 

t- 

X 

* 

IC 

X 

ClC 



cq 

lO 

CM 

id 

CO 

d 

00 

d 

X 

ft 

05 

CO 

CO 

05 

d 

lO 

d 

(35 

d 

ft 

X 

d 

00 

d 

C3 

X 

id 


05 

CM 

(35 

cd 

b-* 

(X 

X 

cd 

X 

d 

X 

rO 

CO 

CO 

03 












ft 












ft 


o 


fe 

ft 

ft 

d 

ft 

ft 

ft 

c5 

ft 

ft 

03 

ft 

ft 

ft 

ft 

03 

ft 

ft 

cS 

ft 

ft 

d 




CO 

CO 

CO 

ft 

CO 

CO 

CO 

ft 

CO 

CO 

a 

ft 

CO 

CO 

CO 

CO 

ft 

CO 

CO 

ft 

CO 

CO 

a 

ft 

d 

CO 

O 

























K 




05 

o 


O 

t- 



X 

CM 


o 

zo 

lO 

CM 


t- 

lO 

ZO 

m 


X 

05 




CM 

o 

o 

CO 

o 

lO 

o 

CO 

CM 

o 

X 


CM 

CM 

(35 

05 

(35 

ZO 

o 

o 

o 

X 

X 

X 


> 

D 


o 

o 

ft 

d 

d 

ft 

CM 

d 

d 

d 


d 

d 

d 

d 

d 

ft 

d 

d 

ft 

d 

d 

4J 

o3 

-H 

-H 

-H 



-H 

-H 


+i 




-H 

+1 

+1 



-H 

+1 

+i 

41 



+1 



o 

lO 

CM 

00 

CM 


t- 


X 


(35 


X 


(35 

b- 


X 

lO 

1^ 

o 

X 

b- 

o 

X 


05 

t- 

00 

CO 

O 

o 


A 

X 

o 

(35 

A 

t- 

X 

1^ 

lO 

A 

I-H 


X 

1-H 

lO 

ZO 

b- 





00 

id 


d 

ft 

cd 

CM 

ft 

cd 

cd 

ft 

id 

CM 

cd 


d 

cd 

d 

cd 

W 



CM 

CM 

CM 

CM 

CM 

CM 

CM 


CM 

CM 

CM 


CM 

CM 

CM 

CM 


CM 

CM 

CM 

CM 

CM 

CM 

CM 




























< 

o 

Exp. 

CO 

10629 

16222 

92341 

3272 

CO 

o 

o 

X 

1670 

206 

lO 

X 

X 

1627 

4954 

4024 

3053 

4130 

1643 

544 

2147 

336 

145 

208 

4063 

2301 

108 

107 

224 




















X 

t- 

X 


CM 


05 

lO 




















X 

X 



O 


O 



. 


















1—1 

1—1 

CM 


CM 


CM 

CM 


CO 


















d 

d 

d 


d 


d 

d 


II 


















-H 

-H 

-H 



ft 


+1 




















lO 

t- 



X 


X 

ZO 




















(.D 

X 

o 


ZO 


(35 

X 




















CM 

CD 

(3D 


q 

O 

X 

q 

I- 



















d 

d 

d 


d 

d 

d 

d 

8cq 














































X 

X 



o 


o 

b- 

<] 

lO 


















b- 

X 



o 


(35 

05 

CM 


















I-H 

I-H 

CM 


CM 


1-H 

1-H 

CM 


















d 

d 



d 


d 

d 

CM 

II 


















-H 

-H 

ft 



ft 


+1 

!/) 




















o 



lO 


b- 

ZO 

Q 



















CD 

lO 



X 


ZO 

X 




















CM 

CD 

O 


q 

o 

CO 

q 

S 



















d 

d 

d 


(d 

d 

cd 

d 




















o 

o 


(35 

05 


X 

CM 




















t- 



lO 

05 


b- 

b- 




















1—1 

1—1 


CM 

1-H 


1-H 

1-H 


CM 


















d 

d 


d 

d 


d 

d 


II 


















-H 

-H 

ft 

+i 


ft 


+1 




















X 



(35 

t- 


ZO 

b- 



















1—1 

t- 


X 

o 


X 

X 




















CM 

o 

o 

1^ 


o 

X 

1-H 




















d 

d 

d 

d 

cd 

d 

cd 

d 







t- 

X 

X 






X 

X 

t- 


t- 

o 

o 

ZO 

lO 

b- 


X 



yCq 







O 








1^ 


o 

ZSJ 


CO 




X 


CO 



rH 





T—1 






’—1 

’—1 

1^ 



1—1 

1—1 

1^ 

1^ 




'“1 

Uh 

<d 





d 

d 

d 






d 

d 

d 


d 

d 

d 

d 

d 

d 


d 

d 

+ 


-H 

-H 

-H 


ft 

-H 

-H 

-H 

-H 



-H 

“H 

-H 

-H 

ft 

-H 

-H 

ft 

-H 




+1 








X 







(35 

X 

o 


X 

X 

b~ 

X 

X 

o 










o 

(35 

(05 







X 

o 



ZSJ 

o 


1^ 

lO 


X 

X 



CM 

o 





o 


o 

o 

o 











o 


ft 


< 


O 

o 

d 

o 

(05 

d 

d 

d 

d 

d 

d 

O 

o 

d 

d 

cd 

d 

o 

d 

d 

cq 

d 

cd 

o 


CO 


o 

CM 




X 

05 

X 

X 


CM 


X 

CM 


o 

o 


o 

(35 

X 

ZO 

X 

lO 




o 

CO 

CO 

CO 


X 

t- 

(X) 

X 

X 

X 

t- 

X 

ZO 

(35 

(35 

lO 

X 

X 

X 


m 

X 

b- 



1—1 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

O 

CJ) 

o 

o 

o 

o 

o 

o 

o 




o 

o 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

flj 






CO 

CO 

CO 

CO 

CO 

CO 

CO 


CO 

CO 

CO 

CO 

CO 

t- 

t- 

t- 

b' 

b- 

b- 

b- 

b~ 

i-. 

CO 



CO 

ft 

CO 

ft 

CO 

ft 

U} 

Q 

CO 

U} 

Q 

CO 

C/J 

Q 

CO 

C/J 

Q 

CO 

C/J 

Q 

CO 

CO 

Q 

CO 

U} 

Q 

CO 

CO 

ft 

C/J 

Q 

CO 

CO 

Q 

CO 

C/J 

Q 

CO 

C/J 

Q 

CO 

C/J 

Q 

CO 

ft 

ft 

X 

ft 

ft 

X 

2 

X 

2 

X 

2 

X 

2 

X 

2 

X 

2 

2 

X 




o 

CO 

CM 








X 










bn 




CD 




CM 


X 

lO 


05 


X 

o 

o 

o 

X 

lO 


(35 

o3 

o 


>5, 

d 

M 

0 

d 

d 

(fl 



o 

o 

CM 

CM 

o 


I> 

X 

05 

lO 

X 


lO 



05 

CD 

o 

o 

O 




CM 



o 

o 

O 

o 

X 


l> 

1-^ 

X 

(35 

o 

CM 

CM 


CM 

X 

05 

(35 

(35 

(35 

(35 

05 

(35 

05 




1^ 

r- 

o 

o 

CM 

CM 

CM 

CM 

CM 

X 

t- 

t- 


X 

05 

05 

(35 

(35 

(35 

05 

05 

05 

05 

Iz; 



o 

o 

o 



1—H 

rH 

1^ 

1^ 





1^ 

1^ 




1^ 

1^ 




f-H 


A 

D^ 


O -2 


< ^ 


w O 


O 


CO X 
CQ 

^ CO 


rJ:3 


^ ” 

Qi :z 

1^ CO 


isophotal radius of the host. 

* Visual inspection found that this SNIa was within the isophotal radius of it’s host. A dust correction was applied. 




































17 


bbbbbbb^^^ 

ododood‘?Oi-H 

bbbbbbbbbb 

dddrHdddddi-H 


ooooi>t-ooa5t>2Soi'^ 

(MCMCQCS|(MC^CM£3 c^(M 

oqoooooPoo 

ddddddd'vi^'^ 

-H^-H^-H-H^HH-H-H 

COCN^OOOOiiO^^Z^OOOi 

--i^i-icsioo^Soco 

0000000*^00 

ddddddd‘9oo 


bbbbbbbbbb 

OOi-HOi-iO^^'-HOO 


bbbbbbbbbb 

C^T—(fOOOCOt^COlOi—li—I 


iO JiO lO »iO 
(M CM <M IM 

o o o o 
d d d d 
-H ^ -H ^ 

O O 1-H 
(M CM fO 

o q q q 
d d d d 


lo io m o lo CM 


0 0 0 0 0 0 3 

-H -H ^ -H -H -H o 

I-I 00 ^ CM (M o r-; 

CO --H CO 1-H o ^ ^ 

q q q q o q 

d d d d d d 


bbbbbbbbbb 

iOOlOOOO'-HCOCO -i-H 

cmcmcmcmcmcmcm3?cm 

bbbbbbbbbb 

OI>O^OCM"^'^OCM 

CMCMCMCOCOCMCM^dcM 


CO CO CO CO 
CM CM CM CM 

o o o o 
d d d d 
-H d -H ^ 

O CM O ^ 
iO O O N- 

q q q q 
d d d d 


CM CO lO i-l 

CM CM CM CM ^ CM 

0 0 0 0 ^0 

d d d d 

-H -H ^ -H -H 

l> <—I t> ^ l> 
O iO uo CO ^ 

q q q q ^ q 

d d d d 9* o 


b'«^bbc^co^^ 

O • CO O • • CM ^ 

33.^33^3^313 

bbbbbbbb 

O'^COOCOCOCM'^ 


oo*^oooo*^'^ooo®^o 

o^oo'^^oo^o 

6°6d°°66°6 

02SlOt>‘5i2CMl>qiO) 

O^O^^^CMO^O 

OHOO'H^OO^O 

d'^dd'^^dd^d 


o o II II 

6 6 II II 3 

q o ^ ^ ^ 

^ CO "O ”0 ^ 

AA § § ^ 

Q g S 


OO COTOC^)^,—IrT 


ttitHCxacA 
0O3 CO COr^r^^ 
^^^WWi—icO"^ 



18 


CD Q 
Q) fl 

2 -J 

03 

H a 

Sh 

tS 


Sh 

o3 

4^ 

CO 




_ 

_ 




_ 

_ 








b 

b 

b 

b 

b 

b 

b 

b 

_ 

b 





iO 

00 

CO 

;o 

CO 

N 

CM 

b- 

b 

O 





ci 


CM 

CM 

CM 


CM 








1 


1 



1 


1 

d 

1 













'—' 






b 

b 

b 

b 

b 

b 

b 

b 

b 

b 





CO 

CO 

O 

CO 

CO 


05 



CO 



CO 

II 

<D 

CO 

1 

CO 

CO 

1 

CO 

CO 

CM 

1 

CM 

CM 

1 

d 

iH 



« 

CO 

h- 

CO 

CO 

CO 

b- 

CO 

b- 


to 




(M 

CN 

CN 

(M 

CN 


(M 

CM 


CM 



rn 

M 

O 

o 

CJ 

CJ 

CJ 

CJ 

CJ 

CJ 

o 

CJ 



CJ 


O 

d 

d 

d 

d 

d 

d 

d 


d 

CM 


hJ 

b 

+1 


+1 

+1 


+1 

+1 

+1 

-H 




§ 


CO 

CO 


b- 

h- 

oo 


b- 


CM 


X 


CX) 


1-^ 

(X) 

tXJ 

lO 

b- 

lO 




o 



o 


CJ 

CJ 

CJ 

CJ 

CJ 

CJ 

o 

CJ 




w 

b 

o 

d 

1 

d 

d 

d 

1 

d 

d 

d 

d 

d 

1 


fi 

cd 


<D 

lO 

UO 

lO 

lO 

lO 

lO 

lO 

lO 

lO 





iz; 

lO 

lO 

lO 

lO 


lO 

lO 

lO 

lO 



> 


CO 



1—1 

1—' 

1 

1—' 


1—1 

1—' 

1 



















S 

V 



V 

y—S 


V 


s 





b 

b 

b 

b 

b 

b 

b 

b 

_ 

b 


Id 



00 


CO 

00 

CJO 

CO 


O 

b 

CO 









d 



05 

d 


0 


hb 

1 

< 

1 

, 1 , 

, 1 , 

1 

, 1 , 

1 


, 1, 


bi; 



b 

b 

b 

b 

b 

b 

b 

b 

b 

b 


o 

lo 



iO 

05 


CM 


N 



00 


d 


ci 


rH 

CM 

CM 

O 

iH 

rH 


o 



II 













6 














0 




>r> 

-O' 

lO 

>r> 

to 

to 

CO 


CO 





CN 

(M 

(M 

CN 

<M 

iM 

CM 


CM 



rn 

M 

O 

CJ' 

CJ 

CJ 

CJ 

CJ 

CJ 

CJ 

o 

CJ 



u 


o 

d 

d 

d 

d 

d 

d 

d 


d 

CO 

— 

J 

b 

+1 


+1 

+1 


+1 

+1 

+1 

-H 



S-i 

§ 


CO 




CO 

05 

CM 

05 


00 




lO 

CO 

-O' 

lO 

lO 

C J 


CM 







o 

CJ 

CJ 

CJ 

CJ 

CJ 

CJ 

CJ 

o 

o 





o 

d 

d 

d 

d 

d 

d 

d 


d 


cd 



1 

1 

1 

1 

1 

1 

1 

1 


1 


T5 














G 














3 


<D 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 





z 

lO 

uo 

lO 

iO 


lO 

lO 

lO 

lO 





CO 

1—1 

I—( 

1—1 

1—1 

I—( 

1—1 

1—1 

1—1 

1—1 

I—( 


d 




























CO 




s 


s 

s 


s 



s 





b 

b 

b 

b 

b 

.—. 

b 

b 

.—. 

b 


i/T 




CO 

1—1 

CM 


b 

00 

CD 

b 

CO 










05 

d 

d 

lO 

d 


cd 






1 


d 



d 











'—' 



'—' 



cd 



b 

b 

b 

b 

b 

b 

b 

b 

b 

b 


O 

tSJ 




iO 

CM 

CO 

CO 


05 

CD 


CO 









i-H 

d 

d 

d 

d 


” 



1 

< 

1 

1 

1 



1 





II 













Si 

3 

Pi 




-O' 





lO 


CM 


o 


rsi 

fN 



fN 


CN 



rM 


rn 

M 

o 

CJ 

CJ 

CJ 

CJ 

O 

CJ 

CJ 

o 

CJ 


cd 

(J 


o 

d 

d 

d 

d 


d 

d 


d 

CO 

hJ 

b 

+1 


+1 

+1 


-H 

+1 

+1 

-H 



bb 

§ 


CO 

>r< 

00 

rsi 



CO 

CO 


00 




CO 

CO 

(M 


CO 


CM 

1—1 


C J 





o 

CJ 

CJ 

CJ 

CJ 

o 

CJ 

o 

o 

CJ 




w 

b 

o 

d 

1 

d 

d 

d 

1 

d 

d 

d 

d 

d 

1 


cd 

’C 


<D 

1> 

h- 

h- 

b- 

h- 

b- 

b- 

b- 

b- 

CO 


> 


z 

lO 

UO 

lO 

iO 


lO 

lO 

lO 

lO 



C5 


CO 

1—1 

1—( 

1—1 

1—1 

1—( 

1—1 

1—1 

1—1 

1—1 

1—1 


4) 














a 































b 

b 

b 

b 

b 

b 

b 

b 

b 

b 





N 


05 

CO 

05 

CM 

CM 

QO 


C5 


cd 

o 



6 

o 

O 

O 

O 

O 

o 

O 

CM 

o 




sd 


sd 

sd 

sd 

sd 

sd 

sd 

sd 




CO 














b 

b 

b 

b 

b 

b 

b 

b 

b 

b 





N 

CM 

05 

CO 

O 

CM 

CM 

QO 


05 


q 



6 

o 

O 

o 

iH 

O 

o 

O 

CM 

o 



(M 













> 

o 


0) 

00 


00 

00 

00 

00 

05 

1—1 

OO 

b- 


-0 










CM 





c 


o 

o 

o 

o 

o 

o 

O 

O 

o 

o 


cd 

CO 


o 

d 

d 

d 

d 

d 

d 

d 

d 

d 

CO 

-0 


b 

+1 


+1 

+1 


+1 

+1 

+1 

+1 






CO 

CO 

CO 

CO 

b- 

-0 

CO 

CO 

00 







C J 

1—1 

c:> 


C J 

C ^ 

1—1 

CC.5 





Ijh 

w 

b 

o 

o 

CJ 

d 

1 

o 

d 

CJ 

d 

o 

d 

CJ 

d 

CJ 

d 

o 

d 

CJ 

d 

CD 

d 

1 


o 

j: 


<D 












S 

0 


Oi 

Oi 

05 

05 

05 

05 

05 

05 

05 

CO 




z 



b- 

b- 

b- 

b- 

b- 

b- 

b- 

b- 




CO 



1—1 





1—1 




fl) 














3 











Pi 



CJ 

S 

0 











d 

CO 








3.1 

2.7 




"Id 



w 




CO 

CO 

II 

1 

II 




o 



Si 




CO 

cc 

W5 






(D 

hJD 












Si 




+1 


>1 

2 



o 





S 




CO 

cd 

cd 

cd 

Cd 




.2 


o3 




T1 

-n 

j 

Sh 


bO 


"cd 

a 

o 

cc 


o 

CO 


uv)= 

ji^ 

p- 

b 

0 

0 

0 

JU 

0 

0 

0 

JIJ 

j 

o 

a 

cd 

<D 

a 

Cd 

0) 

a 

cd 

4P> 

CO 

.B 

rt 

'Bh 

f-, 

0 

f-, 

f-, 

H 




tlH 

fcn 

« 


o 

u 

CJ 

cd 


cc 








iX 

iX 

Uh 


b 


o 


fl 

o 



Cn 

05 

xi 







<; 

2; 

P, 

b, 

w 

w 


CO 


o 

CM 

CO 

d 

td 





ci 

o 

CQ 

Pi 

03 


N 

0 ) 


I 





.2 

2 

i-. 

& 


iS 

03 



b b 


rH 00 

CO 



CD 

o o 

pLi 

o d 

1 -H c/3 

-H-H 

CO t3 

lO CD 

II 

CD iD 

Cd (M 

Pi 

d d 

CO 


O 



o o 

s 1 

m 

n3 

d, 

d d 
-H-H 

X t> 

b 

(73 l> 


d d 

0) 

ID ID 

'Z 

ID ID 

CO 

^ t-H 




c§ q 


d o 


ID O 

CD 

p p 

pin 

d d 

lO C/3 

-H-H 

b 

II 

X t-H 

X l> 


d d 



O 


i-J » 

p p 

§ 1 
m 
cd 
d 

d d 
-H-H 

X C7> 


(73 CD 


1—1 t—1 


d d 

CD 

CO CO 

Iz; 

ID ID 

03 



b b 

1—1 ID 

in 

d d 


CD ID 
^ X 
p p 

pin 

d d 

o c/3 

-H-H 

c<i >= 

II 

i-H (73 
^ (M 

1—1 1—1 


d d 

Pi 


m 

X X 

O 



o o 

§ « 
cd 

d d 
-H-H 

ID CD 




d d 

CD 

b- b- 


ID ID 

03 

1^ 1-H 


b b 

X q 


d o 


^ CD 

X (M 

P p 

pLi 

d d 

C/3 

-H-H 

b 

X X 

<M 

1—1 t—1 

H 

d d 



< 


CO 

o o 

1^ T—H 

0) 

p p 

'm 

m 

cd 

d 

d d 
-H-H 

t- 


CC! I-H 


1—1 1—1 


d d 

CD 

(7 7 


b- b» 

CO 






n3 


X! 


rt 




0 

X! 

X 

D T? 



PlH 

C/3 

W 

b~ X 



19 








